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Abstract 



This paper discusses the construction of high-performance ground-based gamma-ray Cherenkov telescopes with a Davies-Cotton 
reflector. For the design of such telescopes, usually physics constrains the field-of-view, while the photo-sensor size is defined by 
limited options. Including the effect of light-concentrators in front of the photo sensor, it is demonstrated that these constraints 
are enough to mutually constrain all other design parameters. The dependability of the various design parameters naturally arises 
once a relationship between the value of the point-spread functions at the edge of the field-of-view and the pixel field-of-view is 
introduced. To be able to include this constraint into a system of equations, an analytical description for the point-spread function 
of a tessellated Davies-Cotton reflector is derived from Taylor developments and ray-tracing simulations. Including higher order 
terms renders the result precise on the percent level. 

Design curves are provided within the typical phase space of Cherenkov telescopes. The impact of all design parameters on the 
overall design is discussed. Allowing an immediate comparison of several options with identical physics performance allows the 
determination of the most cost efficient solution. Emphasize is given on the possible application of solid light concentrators with 
their typically about two times better concentration allowing the use of small photo sensors such as Geiger-mode avalanche photo 
diodes. This is discussed in more details in the context of possible design options for the Cherenkov Telescope Array. In particular, 
a solution for a 60 mm 2 photo sensor with hollow cone is compared to a 36 mm 2 with solid cone. 

Keywords: TeV Cherenkov astronomy, Davies-Cotton design parameters, photo-sensors, Winston cone 
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The next generation Cherenkov Telescope Array (CTA) will 
be constituted of three types of imaging atmospheric Cherenkov 
telescopes: small-, medium- and large-sized telescopes. The 
array will be dedicated to the observation of the high energy 
gamma-ray sky with unprecedented sensitivity over a broad 
range of energies (0.01 TeV < E y < 100 TeV). This instru- 
ment will enable astronomers and astro-particle physicists to re- 
fine models of gamma-ray sources and underlying non-thermal 
mechanism at work, to study the origin and the composition of 
the cosmic rays up to the knee region, question the nature of the 
dark matter, etc. 

These kind of telescopes detect Cherenkov light emitted 
along developing showers in the atmosphere nearly uniformly 
illuminating the ground over an area of about 50' 000 m 2 . 
Shower maxima of a gamma-ray induced electromagnetic 
shower occur at altitudes comprised between ~7km and 
~12km, for a primary energy ranging from ~100GeV to 
-100 TeV. 

The photon density on the ground depends on the energy of 
the gamma ray, its incidence direction and the distance from the 
resulting shower axis. The displacement of the image centroid 
in the telescope camera is a function of the angular distance of 
the shower core, defined by shower impact distance and alti- 
tude of the shower and is typically about 1.5° for an energy of 
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1 TeV and an impact parameter of 150 m, and 3° to 3.5° for en- 
ergies around 100 TeV with an impact parameter of 300 m. As 
very high energy showers penetrate deeply in the atmosphere 
and generate a large amount of light, they can be observed at a 
relatively large distance from the main light pool already with 
relatively small reflectors translating into the necessity of a de- 
ployment of telescopes with large field-of-view for the explo- 
ration of very high energy showers. In general, the required 
telescope field-of view in order to record contained shower im- 
ages ranges between 3° and 10° depending on the energy range 
of interest. 

A promising design for wide field observations is the Davies- 
Cotton telescope |[T7l . Its advantage is a point-spread function 
enabling a larger field-of-view than a parabolic design. The 
Schwarzschild-Coude design (e.g. ifTBI ) is not considered fur- 
ther, as it presents technical and cost challenges compared to 
a conservative CTA proposal using a prime optics design (see 
also the comments in section [ZT) . 

It is demonstrated that once the field-of-view of the camera 
pixels and of the whole instrument together with the photo- 
sensor technology is fixed by means of physics arguments, only 
one parameter is left free. This parameter can be fixed as well, 
relating the worst optical resolution in the field-of-view, i.e. at 
the edge of the field-of-view, with the size of a pixel. Having a 
reasonable cost model at hand, even the most cost-efficient sin- 
gle telescope or telescope design operating in an array can be 
derived. 
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Knowing the anticipated angular size of a pixel from physics 
constraints, a requirement on the point-spread function of the 
telescope can be determined. If the point-spread function is 
identified by the root-mean-square of the light distribution, re- 
quiring a pixel diameter four times the root-mean-square en- 
sures that most of the light from a point source at infinity is 
concentrated in a single pixel at the edge of the field-of-viewQ 
If simulations show that such a containment is not necessary 
in terms of angular resolution of the shower origin or back- 
ground suppression, also a smaller value like twice the root- 
mean-square can be considered, c.f. (H. If this requirement is 
not optimized or not met, either the physics outcome is wors- 
ened or the camera has more pixels than necessary and will not 
be cost-efficient. 

The CTA array layout will consist of large size telescopes 
(LST, primary mirror diameter ~24 m) in the center, and suc- 
cessively surrounded with an increasing number of medium- 
sized (MST, ~12m) and small-sized (SST, ~6m) telescopes in 
order to instrument a ground surface area comprised between 
4 km 2 and 10 km 2 . Telescope spacings, sizes and field-of-views 
reflect the energy range to be explored by a certain type of tele- 
scope. The LST sub-array will be primarily focusing on the 
observation of the high energy gamma-ray sky with great pre- 
cision below about 100 GeV, the inter-telescope spacing will be 
small, about 60 m and the single telescope field-of-view limited 
to about 3° -4°. The SST sub-array will conversely be opti- 
mized for multi-TeV observations, up to around lOOTeV al- 
lowing for large spacing of up to 300 m -500 m with relatively 
modest reflector size. The camera pixel field-of-view of these 
different telescope types will span between about 0.08° -0.1° 
(LST) up to 0.2° -0.3° (SST), linked to the different shower 
intensity with its intrinsic fluctuations and the concurrent ne- 
cessity of keeping the number of camera pixels reasonably low. 

Consequently, this study will focus on reflector diameters in 
the range of a meter to about 30 m and on pixel field-of-views in 
the range between ~0.08° and ~0.3°. and is extended to off-axis 
angles of the incoming light up to more than 5°. 

In the following, first a description is derived for the rela- 
tions between the existing design parameters of a Cherenkov 
telescopes, e.g. focal length and reflector diameter. Then, by 
including a semi-analytical treatment for the optical quality of 
a generalized Davies-Cotton reflector, this description becomes 
applicable for the design of real telescopes. At the end, the 
influence of the variation of the parameters on the optimized 
design is discussed. 

2. Optimization 

The light collection is one of the most important parame- 
ter of a Cherenkov telescope and can be improved by, e.g., an 
increase of the photo-detection efficiency of the photo-sensors 
or a rescaling of the system, i.e., a corresponding rescaling of 
the reflective area and photo-sensor size. As current technol- 
ogy (photo-multiplier tubes, hybrid photo-detectors or silicon 
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photo-multipliers) only offer a limited choice of photo-sensor 
sizes, the Cherenkov telescope design parameter phase space is 
reduced. Given a particular photo-sensor type, the addition of 
a light concentrator in front of the photo-sensor is the only way 
to increase the light collection efficiency. 

In the following, it will be shown that the relationship be- 
tween the characteristic parameters of an optimal telescope de- 
sign (pixel field-of-view and linear size, reflector diameter and 
focal distance) is fully constrained, once the technological re- 
quirement (photo-sensor size and the light-guide material) and 
the physics requirements (the pixel angular size and the field- 
of-view) are frozen. This result is obtained by imposing restric- 
tions on the value of the point-spread function at the edge of the 
field-of-view in order to keep adequate image quality and thus 
analysis potential. 

2.1. Light concentrators 

The theorem of Liouville states that the maximum concentra- 
tion theoretically achievable is defined by maintaining the phase 
space, i.e. the product of solid angle, defined by the incoming 
light rays directions, light ray momentum squared and the sur- 
face area crossed by the light rays. The theorem of Liouville 
is applicable to the case of a Winston cone with entrance area 
A placed in front of a photo-detector of area A', corresponding 
to the exit area of the cone. A is defined provided A', the solid 
angle Q defined by the light rays entering the cone and the solid 
angle Q.' defined by the light rays leaving the cone. 

Winston has shown in |6l that the maximum concentration 
factor for a rotationally symmetric light concentrator is 
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where n and n' denote the refractive index of the media in front 
of and inside the optical system with n w 1 in air. <p corresponds 
to the maximum angle at which a light ray enters the system 
related to the solid angle Q. 

If the system is not axisymmetric or the angular acceptance 
of the photo detector is smaller than £2' = 2n (as assumed in 
Eq. [TJ, C max has to be adapted accordingly. 

Besides the increase of the light collection area, the use of 
cones enables a partial screening of the night-sky light pollu- 
tion corresponding to (f> being larger than the angle of light rays 
coming from the edge of the reflector. 

Simulations [7 1 and recent concentration efficiency measure- 
ments [ 8 1 of solid cones («' w 1 .4) designed for the FACT cam- 
era |0 demonstrated that their shape is nearly ideal, that the 
concentration factor reaches a value close to C max and the ge- 
ometric loss is only of the order of a few percent excluding 
absorption loss^j In the case of hollow cones, Fresnel reflection 
losses have to be considered at the surface of the photo-sensor. 
If the camera is sealed with a protective window, which is usu- 
ally the case, also losses at the window surface need to be taken 
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into account. By choosing a material for the cones and the pro- 
tective window with a similar refractive index than the material 
of the photo-sensor light entrance, these losses can be omitted. 
Combined with the almost perfect reflectivity of solid cones due 
to total reflection (limited only by the surface roughness), solid 
cone usually outperform hollow cone. 

The concentration factor achieved with the Winston cones is 
fundamentally similar to the size reduction of the focal plane in 
a Schwarzschild-Coude design and linked to the conservation 
of the space-momentum phase space according to the Liouville 
theorem: the conversion of the spatial into momentum phase 
space, by means of Winston cones or secondary optics. While 
cones reduce the acceptance of the incoming light rays from 
a large area at the entrance to a large angular acceptance and 
small area at the cone exit, the secondary mirror optics leads to 
a similar spatial compression and angular widening of the light 
rays at the photo-sensor and thus a reduced plate-scale. While 
cones are non-imaging devices, the secondary optics is imag- 
ing. Hence, in the Schwarzschild-Coude design, it is possible to 
attain excellent optical resolution, in terms of Cherenkov tele- 
scope requirements, with a field-of-view as large as 15°. Both 
designs enable compression of the photo sensitive area by fac- 
tors larger than ten w.r.t. their primary optics design. 

2.2. Connection to the optical system 

In the case of a Cherenkov telescope, the light entering the 
cone comes from a reflector visible under a maximum angle 
<p ^ The opening angle of the light at the entry of the cone is 
therefore well defined by the properties of the optical system, 
i.e. by the diameter of the reflector D and the focal length F, 
f = FID. 

tan0=— (2) 
Combining this with Eq.[TJyields 
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For instance, taking the FACT values, n' = 1.4 and / = 1.4, 
we obtain the theoretically maximal achievable concentration 
factor C max = 17.3, i.e. the linear size of the entrance area can 
ideally be larger than four times the linear photo-sensor size. 

2.3. The optical system 

In addition, the optical system defines the zoom factor or 
plate-scale, i.e. the field-of-view corresponding to a physical 
area in the focal plane. The correspondence between the angu- 
lar size § and the linear size 6 on the focal plane is 
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or in the limit of small 



6~$F. 
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Cameras in Cherenkov telescopes are pixelized due to the 
use of photo-detectors. To increase the light collection effi- 
ciency further, and to maintain symmetry, these pixels are usu- 
ally aligned on a hexagonal grid, i.e. in closed package geom- 
etry. In recent years, MAGIC has exploited the photon arrival 
time extracted from the measured pulse and demonstrated sig- 
nificant improvements in the sensitivity ifTUl . 

The technique, taking into account the change of the arrival 
time between neighboring pixels, performs best, if all neigh- 
bors are at an identical distance from the central pixel. Conse- 
quently, the ideal shape of a pixel is hexagonal. 

The distance 6 on the camera surface is the distance between 
two parallel sides of a hexagon, its area is 



A = 
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Combining with Eq.|3] the plate-scale formula Eq.|4]and Eq. |6j 
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is obtained, which translates the close relationship between the 
pixel field-of-view §, the focal length F and the reflector diam- 
eter D, once the technological parameters fixed: photo detector 
size A' and light concentrator material n' . 

Defining a constant related to these properties 



k(n,n',A') = ^.(^) 



and rewriting Eq.[7](tani? » & for § <sl°) as 

1 



F = 
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it is immediately apparent that the focal length F of the sys- 
tem is a direct consequence of the pixel field-of-view and the 
reflector diameter, if the properties of the photon detector and 
the material of the cones are known. For typical Cherenkov 
telescopes, F/D is between unity and two. Below unity the res- 
olution becomes too coarse and above two, not only the number 
of pixels and hence the price of a camera becomes too high, but 
also the camera holding structure becomes mechanically com- 
plex and hence disproportionally expensive. This constraint on 
F/D applied to Eq.|9]yields 



4.25 < k& 2 D 2 < 5. 



(10) 



> defined from the focal plane center 



Precisely the choice of F and D defines the optical quality 
of a mirror system. At the same time, the size of a single pixel 
defines a natural constraint on the optical quality of a system, 
i.e., F/D should be chosen such that the point-spread function 
at the edge of the camera is within a limit well defined by the 
pixel's field-of-view. 

2.4. Optical quality 

The light collection area is important for a Cherenkov tele- 
scope and typical reflector sizes range from a few to ten or 
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twenty meters. However, with the current technology, it is not 
possible to produce large mirrors with the requested quality at 
a reasonable cost. Furthermore, optical systems compiled from 
a single mirror suffer large aberration effects at large off-axis 
angles, while a wide field-of-view is necessary for the obser- 
vation of multi-TeV showers up to large impact parameters, as 
well as for extended sources. Therefore, segmented mirrors are 
in use. The layout providing the best optical quality for seg- 
ments of identical focal length is the so-called Davies-Cotton 
layout [ 15 1, where the single spherical mirrors are located on a 
sphere with radius F and focused to a point at 2F. 

The relevant quantity which influences the on-axis and off- 
axis optical quality is the focal ratio F/D. The optical quality 
improves with larger values. This scale invariance statement is 
true only as long as the optical quality of a single mirror can be 
neglected against the optical quality of the whole system, which 
is generally the case at the edge of the camera. 

To be able to constrain the optical point-spread function, a re- 
lation between the tessellation, the focal ratio and the resulting 
point-spread function is needed for a given maximum inclina- 
tion angle of the light, i.e. at the edge of the field-of-view of the 
camera. 

In Appendix A.l a formalism describing the point-spread 
function of an ideal Davies-Cotton reflector, i.e. a reflector 
with infinite tessellation is presented. The point-spread func- 
tion is described by the root-mean-square of the light distri- 
bution. With the help of ray-tracing simulation, a reasonably 
good description of a real tessellated Davies-Cotton reflector is 
derived from this analytical approach in Appendix A. 2 In- 



cluding the correction factor which describes the deviation of 
the analytical approach from the simulations, a good descrip- 
tion of the point-spread function is obtained. It is shown that 
the point-spread function P of an ideal Davies-Cotton can be 
expanded into a polynomial in a' ■ f~K Given the incident angle 
a <K 1 of the incoming ray and / in the range between one and 
two, the polynomial is hypothesized into the more simpler form 



P(J,a,N) = c (a,N)-r 



(11) 



with coefficients cq and c\ and N the tessellation number as 



described in Appendix A. 2 This parametrization is found to 
match the ray-tracing simulation without loss of precision. The 
coefficient co can directly be deduced as the result of Eqs. A. 1 1 
at / = 1 and c\ is derived by a fit. An example for the coeffi- 
cients co and c\ for selected a is shown in Fig.[T] 

2.5. Result 

As discussed in the introduction, it is required that the point- 
spread function is small compared to the pixel field-of-view at 
the edge of the field-of-view, so that the light of a point source is 
well contained in one pixel. Defining the ratio r between both, 
this requirement can be expressed as 



§ = r-P. 
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Combined with Eq.[TT] the focal ratio / can now be expressed 
as 



/ = 



A 



r ■ c 



(13) 



Including this in Eq.[9]yields 
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with a correction factor g defined as 
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The absolute focal length F can now be calculated using Eq.|9] 
To deduce the effective reflective area from Eq. [14J the 
shadow of the camera on the reflector has to be taken into ac- 
count. To calculate the fraction of the reflector shadowed, the 
ratio of their areas is calculated. Conversion of a from an angle 
to a length yields approximately R = F ■ a for small values of a 
(«10°). 



Expressing the focal length by Eq. 13 the fraction of the cam- 
era shadow on the reflector is derived. 
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If a real camera housing is significantly larger than the photo 
sensitive area itself, a correction factor should be included. 

Now the effective light collection area of the optical system 
can be deduced as 



A eff = 



Kg 
k 



1 - 
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(17) 



If real setups should be compared, like e.g. Davies-Cotton 
and Schwarzschild-Coude, also other sources of light-losses 
must be included, such as geometrical efficiency of the cones 
(light-loss at the edge of the mirror), total mirror reflectivity, 
cone transmission or reflection losses, or photo detection effi- 
ciency. 

2.6. Discussion 



The relation given in Eq. 14 includes several parameters 
which are subject to change. For simplicity, a standard setup 
has been defined to which altered setups are compared. 

Silicon photo-detectors are a recent and very promising tech- 
nology. Therefore, a silicon photo-detector with a sensitive 
area of 36 mm 2 is chosen as a benchmark device. Such devices 
are already commercially available with acceptable properties. 
Although their sensitive area is still rather small compared to 
photo-multipliers, by increasing their light-collection with solid 
light concentrators, their light-collection area becomes reason- 
ably large. Such light-concentrators still maintain a reasonable 
weight and length in term of absorption. Typical Plexiglas ma- 
terials have a refractive indices of the order of n'=lA and are 
used hereafter as a reference. 

As the light-collection area of a telescopes scales directly 
with the photo-sensitive area, the most obvious use of small 
photo sensors is a small telescope sensitive mostly to high en- 
ergetic showers. At high energies, the collection area of a tele- 
scope array is of prime importance due to rapidly decreasing 
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Figure 1 : Coefficients derived from the work in the Appendix to convert the focal ratio / into point-spread function as defined in Eq. |ll| Sagittal component (left) 
and the tangential component (right). The dashed line denotes on-axis rays a=0°, the solid lines (from the dashed line outwards) rays at 2°, 4° and 6° off-axis. 



fluxes. Due to the bright light-pool of high energy showers, 
telescopes with relatively small reflectors can be operated with 
large spacing of, e.g., 400m or more, c.f. [2|. Since such spac- 
ings demand a large camera field-of-view, a field-of-view of 9° 
diameter is chosen as a reference. 

A typical reflector for a Cherenkov telescope with Davies- 
Cotton layout enables the manufacturing of a primary reflector 
tessellated into spherical mirror of identical focal lengths. From 
the scaling with the tessellation number as derived in |A.l 1| it 
can be concluded that a layout with only three mirrors on the 
diagonal (N=3) has still a significantly worse optical quality 
than a reflector with five mirrors on the diagonal. Although 
the point-spread function at the center of the camera is clearly 
dominated by the mirror size, the relative influence almost van- 
ishes at higher off-axis angles. Since the solution with A^= 3 
still shows a degradation of more than 10% compared to the 
solution with A^= 5 even at the highest simulated off-axis an- 
gles, it is discarded. On the other hand, a further increase of the 
tessellation number (individual mirror size over primary reflec- 
tor diameter) does not significantly improve the optical qual- 
ity. Consequently, choosing A^= 5 is a good compromise and 
already close to the optimum achievable. Comparable results 
were obtained in [4| although using a third order approxima- 
tion overestimating the optical quality. 



It must be noted that the simulation does not take the point- 
spread function of the individual mirrors nor any possible mis- 
alignment into account which must be added quadratically to 
the result. However, for the solutions discussed here this can 
be neglected, c.f. 0. In general, alignment errors can be kept 
minimal and individual mirrors can be machined with a point- 
spread function small compared to the point-spread function at 
the edge of the camera. 

On average, all Davies-Cotton designs with a reasonable 
F/D have a root-mean-square of the light distribution in the 
tangential direction about two times larger than in the sagittal 
direction. 

Ideally, the sagittal root-mean-square at the edge of the cam- 
era should fit a fourth of the pixel's field-of-view. This ensures 
that in the sagittal direction 95% of the light is contained within 
one pixel diameter and roughly 68% in the tangential direction. 
However, since the point-spread function is not Gaussian and 
has long tails in tangential direction exact numbers for the light 
content might slightly differ. 

For convenience, all following plots show dots for F/D = {1, 
1.25, 1.5, 1.75,2}. 

Fig. [2] shows the reflective area versus the pixel's field-of- 
view for comparison in the standard case with and without 
shadowing for different camera field-of- views. Since the effect 




Figure 2: Reflective area (left) and the effective reflective area (right), i.e. including the camera shadow, for the standard setup as discussed in section |2~6] The 
result is shown in the range between F/D = [1, 2]. The dots denote intermediate results at F/D = {1.25, 1.5, 1.75}. For convenience the result is also shown for a 
camera field-of-view of 5° and 12°. The legend gives the corresponding pixel field-of-view i? with its relation to the optical point-spread function P at the edge of 
the camera, the photo sensitive area A' of the photo sensor, the refractive index n of the light collector and the tessellation of the reflector expressed as the number 
N of individual mirrors on its diagonal. 



is comparably small and the mirror diameter is more expres- 
sive, in the following all plots show the mirror diameter rather 
than the reflective surface in the non-obstructed case. 

The effects of changing different input parameters w.r.t. to 
the previously described benchmark configuration are shown in 
Fig. [3] and discussed below. 

Changing the camera field-of-view (Fig. [3] top plot). Changing 
the camera's field-of-view basically shifts the valid range along 
the line, i.e. the range corresponding to F/D = [1.0, 2.0]. That 
means that it is possible to build telescopes identical in optical 
quality, pixel's field-of-view and mirror diameter, but different 
field-of-view resulting simply in a different focal length of the 
system. In short: Changing the field-of-view only changes the 
focal length. 

Changing the optical quality (Fig. [3] middle plots). A change 
in the requirement on the optical quality r directly influences 
F/D, and therefore also shifts the range of reasonable F/D al- 
most linearly in § (left plot). Changing the tessellation (right 
plot) is like changing the requirement on the optical quality. 
While the difference in optical quality between a Davies-Cotton 
layout with three mirrors on the diagonal and five mirrors is still 
significant, all other layouts give identical results within a few 
percent. In short: Any tessellation number >5 gives similar 
results. Changing the requirement on the optical quality only 
changes the focal length. 

Changing the photo sensitive area (Fig. [i] bottom left plot). 
Since the constant k is directly proportional to the size A' of the 
photon detector, the mirror area is directly proportional to the 
size of the photo sensor. If the size of the photon sensor is lim- 
ited, a simple way to increase the field-of-view of a single pixel 
is to sum the signal of several photon counters to a single signal. 
To maintain a hexagonal, i.e. most symmetric layout, summing 
the signal of three, four or seven photon sensors seems appro- 
priate. In short: Assuming an optimized light-concentrator, the 
photo sensor's physical size defines the scale of the system. 



Changing the light concentrator (Fig. [3] bottom right plot). 
Another way to increase the reflective area is an increase of 
the refractive index of the light concentrator entering quadrat- 
ically. Using solid cones made from a Plexiglas material with 
a typical refractive index in the order of 1 .4 allows to increase 
the achievable reflective area by a factor of two compared to 
hollow cones. Since the length of a typical light concentrator 
for an exit of 1 mm diameter is in the order of 3 mm -4 mm, 
weight and light-attenuation, which is dependent on the length 
of the material crossed, will define a natural limit on the sen- 
sor size for which a solid cone is still efficient. For comparison 
reasons not only solid (n'=1.4) cones but also intentionally less 
efficient hollow cones (n=1.0) are shown. Non-optimum hol- 
low cones are typically used in current Cherenkov telescopes, 
in which the sensitive area of standard photo-detectors (PMTs) 
is not a limiting factor. In short: Increasing the refractive index, 
quadratically increases the reflective area of the system. 

Another interesting aspect for the final performance of a tele- 
scope is the collection of background photons from the diffuse 



night-sky background. Here, Eq. 17 leads to an interesting con- 
clusion. Since the rate of the night-sky background photons per 
channel scales with the effective reflective area and the solid an- 
gle corresponding to the field-of-view of the pixels, the night- 
sky background rate R is proportional to 
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(19) 



* r g-h 

For the range of F/D = [1,2], Eq[T3"1 and Eq. [15] yield a cor- 
rection factor g between 1.0625 and 1.25. With them, Eq. \19\ 
can be transformed into R oc c/k, with Eq. |8]into R oc c'A'nK 
It is immediately apparent that the night-sky background rate 
scales with the physical entry area of the pixel. Assuming only 
reasonable camera field-of-views between 3° and 13° diameter, 
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Figure 3: Reflector diameter versus pixel's field-of-view for different setups. The results are shown in the range between F/D= [1,2]. The dots denote intermediate 
results at F/D = {1.25, 1.5, 1.75}. For convenience the result is also shown for a camera field-of-view of 5° and 12°. 
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the coefficient c' = 2/ V3c is between 9.6 and 14.2. This can 
be interpreted such that the night-sky background rate per pixel 
can be considered constant within +10% in the first order along 
the lines of an optimized telescope. The dependence of d on 
F/D and the camera field-of-view is shown in Fig. [4] 




Focal ratio F/D 

Figure 4: The color coded coefficient c' describing the scale of the night-sky- 
background rate versus focal ratio F/D and camera field-of-view (diameter) for 
fixed cone material and photo detector size. 

Quantitatively the night-sky background rate tnsb is given by 

>"NSB(rpixei, A milT01 -, x(v), /(v)) = 7!r 2 A mirror j dv ^-^(v) 

where x(v) an d f( v ) are the photo-sensor's photo-detection ef- 
ficiency and the night-sky background intensity, respectively. 
For silicon photo detectors (MPPC H21), as used in the FACT 
camera, and the night-sky background at La Palma ff3l . tnsb is 
about 150 MHz (much larger than the device dark count rate) 
given a reflective area of 10 m 2 and a pixel field-of-view of 
0.2°. If a cutoff in the photo-detection efficiency is introduced 
at 650 nm (PMT-like behavior), this can be further reduced. In 
general, the night-sky background rate is not a main problem 
in Cherenkov astronomy, as the combined trigger requirements 
of signal among nearest pixel neighbors and within a short time 
lead to its very efficient suppression. 

General considerations. Existing telescopes are usually under- 
designed, i.e. the photo-detectors are larger than necessary or 
the light-concentrators do not reach the maximum possible con- 
centration. 

On the contrary, currently so-called Silicon Photo-multipliers 
have proven their potential in Cherenkov telescopes ifTTll . These 
silicon based photo-detectors usually have a very limited area, 
but used in an optimized setup, their effective physical light 
collection area, i.e. entry of the light concentrator, can be much 
larger. If prices of photo-sensors are compared, this has to be 
taken into account. Not the price per mm 2 physical sensitive 
area, but the price per cone entrance area or field-of-view, has 
to be considered. If cheap enough, the signal of several photo- 
detectors, equipped with individual light concentrators, could 
even be summed. 



For a comparison, in terms of effective reflective area, the 
transmission losses of solid cones and their gain from avoiding 
Fresnel reflection has to be taken into account, as well as the 
reflection losses of hollow cones. 

A note on timing. For an ideal Davies-Cotton reflector, the ar- 
rival time distribution of an instantaneous parallel beam flash 
is practically flat. More precisely, it is linearly decreasing, but 
looks flat on the small interval. The number of photons in the 
arrival time interval [T, T + 6T] is N(T, 5T) cc T. Its width 5T is 
given by D/c ■ f DC (Q.5, 0), where c is the speed of light. The in- 
terval is of the order 1.1 ns for a 4m class reflector (F/D ~ 1.5, 
up to slightly less than 4.5 ns for a 12 m reflector consider- 
ing F/D ~ 1.2. This short time spread is not a problem for 
the observation of showers with a small size telescope as it is 
still small compared to the Cherenkov light flash duration. For 
medium and large size telescopes, a slightly different mirror ar- 
rangement should be chosen if time spread matters. By a mirror 
arrangement, intermediate between a spherical (Davies-Cotton) 
and a parabolic design, the time spread can considerably be 
improved, maintaining the point-spread function almost com- 
pletely. While the point-spread function is dominated by the 
majority of the mirrors, i.e. outermost mirrors, the time spread 
is dominated by the ones with the largest DeltaT mirrors, i.e. 
innermost mirrors. Consequently, moving the innermost mir- 
rors closer to a parabola immediately improves the time-spread 
while the effect on the point-spread function is rather limited. 
Ideally, mirrors on a parabola with adapted focal lengths are 
used, but might be a cost issue. With adapted focal lengths, all 
mirrors are placed at correct focal distance, so that, a similar 
point-spread function than for the Davies-Cotton arrangement 
can be expected. 

Remarks about CTA. Recent results of FACT ifTTI show that a 
reflector in the order of 3.5 m diameter can give already reason- 
able physics results with current analysis and detector technol- 
ogy. Therefore, a 4m diameter reflector for SST is assumed. 
For physics reason, the field-of-view is supposed to be between 
9° and 12° (leading to a reasonable F/D between 1.5 and 1.8 
assuming an optical quality of 4 and a tessellation number of 
5). Requiring a pixel field-of-view in the order of 0.26°, pos- 
sible solutions could be solid cones with a 36 mm 2 G-APD or 
hollow cones with a 60 mm 2 G-APD, see also Fig. [5] 

The manufacturing of 60 mm 2 G-APDs is under discussion 
with Hamamatsu. A rough estimate shows that a solid cone 
for such a device would be about three times longer than for a 
9 mm 2 as used in FACT. Considering the transmission loss of 
10% in the FACT cones [8|, which is a very conservative esti- 
mate, such cones would have a loss in the order of 35%. Since 
solid cones avoid the loss from Fresnel reflection at the sealing 
surface and the G-APD surface, the real light loss would only 
be around 27% assuming that the hollow cone has a reflectivity 
of 100% which in reality is not true. 

Keeping the pixel field-of-view constant, the gain in reflec- 
tive area corresponds to the refractive index of the cone material 
squared. In the case of a refractive index of typical Poly(methyl 
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Figure 5: Reflector diameter versus pixel's field-of-view for selected setups. The shown setups were chosen such that a small-size, medium-size and large-size 
telescope is considered. The horizontal black lines are to guide the eye to possible setups like D=3.5 m or D=5 m, D=12m and D=24m. The displayed range is 
between F/D = [l,2]. The dots denote intermediate results at F/£> = {1.25, 1.5, 1.75}. The Greek letter Z denotes the number of photo sensors compiled into a 
single pixel. 



methacrylate) PMMA of 1.4 this is a gain of «100% reflective 
area, which outperforms the transmission loss significantly. 

Assuming that the manufacturing of a 36 mm 2 G-APD would 
be as easy as of a 60 mm 2 G-APD, one can compare a solution 
with a 36 mm 2 G-APD and a solid cone and a 60 mm 2 hollow 
cone (assuming perfect reflectivity). In this case, the transmis- 
sion loss of the solid cone is around 13% compared to 8% Fres- 
nel loss for the hollow cone. On the other hand, the solution 
with the hollow cone yields a 15% smaller reflective area (same 
pixel field-of-view) or 27% more pixels (same reflective area). 

Assuming further that the price of the camera scales with the 
price of each channels, a reduction of the number of channels 
by almost 30% reduces the costs for the camera significantly. 
Since the costs are also dominated by the price for the photo- 
detectors, and the price of G-APDs, in the first order, scales 
with the sensitive area, it can be estimated that the price for the 
36 mm 2 G-APDs would be almost a factor of two lower than 
for the larger ones. 

In Figure [5]possibly solutions for MST and LST designs are 
shown using G-APDs and solid cones. On both cases it is con- 
venient to sum at least three, or even seven, pixels into one read- 
out channel to keep the ratio F/D low for construction reasons. 



3. Conclusion 

The Davies-Cotton design with its simplicity as compared to 
non validated dual optic systems is assuredly a good option for 
a wide field-of-view, up to 10° -12°, high energy Cherenkov 
telescope. 

With this study, it is possible to scan a wide phase space of 
the design of Cherenkov telescopes or telescope arrays. This 
was achieved by a description of the optical performance of 
Davies-Cotton reflectors and introduction of the effect of light- 
concentrators. In particular, this study provides an analytical 
description of the optical performance of a tessellated Davies- 
Cotton reflector precise enough to enable performance studies 
without the need for dedicated simulations. 

By including the effect of the light-collector into the system 
of equations, the available phase space of design parameters 
is reduced to a single parameter, once the photon detector has 
been chosen and either the pixel field-of-view or the camera 
field-of-view has been fixed by physics constraints. While the 
choice of photo sensor is usually defined by the availability on 
the market, constraints on the camera field-of-view are a result 
of the physics targets. 
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Appendix A. Parametrizing a Davies-Cotton reflector 
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Figure 6: Example design overview for a small size telescope calculated for 
an ideal Davies-Cotton reflector, i.e. N ~ 9, a pixel size four times the point- 
spread function and solid cones with n = 1.4. For a given F/D and camera 
field-of-view, the corresponding reflector diameter and the pixel field-of-view 
can be read. 

If these two parameters are fixed, the whole available phase 
space of possible solution can now be scanned by changing a 
single input parameter. It can, for example, be convenient to 
scan a reasonable range of the focal ratio F/D and derive all 
other parameters accordingly. From the result, the most cost 
efficient solution, or the one performing best in sense of physics 
targets can be chosen. 

For the Cherenkov telescope array, several design options 
were presented. It could be shown that for the small size tele- 
scope, considering a camera field-of-view of 9° to 12°, a four 
meter reflector is enough if 36 mm 2 senors are topped with solid 
cones to achieve a pixel field-of-view in the order or 0.25° to 
0.3° at reasonable F/D. An alternative solution are hollow 
cones with correspondingly larger sensor area, which is disfa- 
vored because of the costs dominated by the sensor. An ex- 
ample plot which easily allows to determine reasonable options 
from the available phase space is shown in Fig. [6] The reflector 
diameter can easily be re-scaled linearly with the photo sensor 
size and the refractive index of the cone material. 

For the medium size and large size telescopes, the most 
reasonable solution using small sensors would be the summa- 
tion of three and seven, respectively. Equipped with differ- 
ent sum-stages, these modules could be applied in any tele- 
scopes. Larger silicon based sensors, expected soon on the 
marked, would allow a single-channel/single-sensor solution. 
Using several small sensors in one channel has the advantage 
that the application of solid cones is possible in terms of weight 
and transmission and costs for photo sensors can be kept low 
due to their at least two times higher compression ratio. 



Appendix A.l. Ideal Davies-Cotton reflector 

The Davies-Cotton design |[T5l is known to be promising 
for wide field prime-focus telescopes and was studied earlier 
analytically |[T6ll and through simulations [IT]. However, the 
parametrizations are moderately accurate and non existing for 
tessellated reflectors. 

Here, parametrizations are provided, accurate at the percent 
level up to 12° field-of-view, for the ideal (non-constructable) 
Davies-Cotton telescope and accurate to a few percent for a re- 
alistic Davies-Cotton telescope with arbitrary tessellation of the 
reflector. 

Prime-focus telescope design. The major issue of the design 
of a telescope is the reflector and its optical performance. Since 
design parameters like the field-of-view of a single pixel or the 
field-of-view of the whole camera are closely related to the re- 
flectors optical performance, it is important to understand the 
relation between the reflector design and its performance. Un- 
fortunately, neither spherical nor parabolic mirrors can provide 
both, good optical point spread function for on-axis and in- 
clined rays, at the same time, because the distance between any 
point on the mirror surface to the focal point does not match the 
local focal length defined by the local radius of curvature. Fur- 
thermore, in the case of a spherical mirror, also the shape of the 
mirror surface is not ideal compared to a parabolic mirror. The 
parabolic shape ensures that parallel rays from infinity are well 
focused into a single point (due to the definition of a parabolic 
surface) while in the spherical case this is not the case. That 
means that in both cases rays hitting the mirror far off its center 
have their focal point not at the focal plane. In the case of a 
spherical mirror they also miss the focal point (aberration). 

Consequently, the ideal mirror would be a combination of 
two properties: A mirror surface which is shaped such that it 
has the right focal distance at any point, but at the same time 
any point is correctly oriented, so that focal distance and direc- 
tion are correct. Since local normal vector and local curvature 
cannot be disentangled such a mirror can only be a theoretical 
construction. Tessellating the reflector into individual mirrors, 
this behavior can be approximated, as shown by Davies and 
Cotton, if the reflector is built from several spherical mirrors 
which are placed on a sphere around the focal point. In this 
case, the reflector can have the correct focal distance locally 
and, at the same time, the mirror elements can be oriented such 
that they correctly focus to the focal point. Apart from an im- 
proved optical performance for inclined rays, the production of 
several small and identical mirrors is also much more cost effi- 
cient than the production of a single large mirror. 

Since any optical system can always be linearly scaled, in the 
following a scale factor is chosen such that the reflector diame- 
ter corresponds to unity, which is identical to defining f — F/D 
with F being the focal length and D the diameter of the mirror. 

Spherical reflector. The spherical mirror has its focal point at 
half its radius of curvature /. Its surface is given by z — 
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f S pb(x,y) and its normal vector by n sp h(x, y): 

f spb (x,y) = If- ^(2f) 2 ~(x 2 +y 2 ) (A.l) 
» sp h(-*> y) = (x/(2f - / sph ), y/(2/ - / sph ), - 1) (A.2) 

Ideal Davies-Cotton reflector. The ideal Davies-Cotton reflec- 
tor has a non constructable surface. Its shape z = foc(x,y) is 
spherical with radius of curvature /, but its normal vectors are 
defined to intercept at location F = (0, 0, 2f). Formally, the sur- 
face equation foc(x,y) and the surface normal vectors fiDc(x,y) 
are 



foc(x,y) 
n D c(x,y) 



f- ylf 2 -(x 2 +y 2 ) 
WDc(x,y)-z) 



(A3) 
(A.4) 



y 



2/ - foe (x,y) 2f- foe (x, y) 



Practically, this equations describes infinitely small mirror 
elements placed on a sphere, oriented accordingly. 

Taylor development. To have the root-mean-square of the pro- 
jection of reflected rays on the focal plane along x and y coin- 
ciding with tangential and sagittal resolutions, a rotation around 
z is performed without loss of generality. 

An incoming ray with vector v = (0, sin</3, cos0) will 
therefore be reflected on the surface in the direction v r - 
v - 2(v • n)njn 2 and intercept the (non curved) focal plane at 
r = -y/x 2 + y 2 , generally yielding 



(X, Y,Z) = \x+ ^(/ - f(x,y)), y + ^(/ - f{x,y)\ f 



(A.5) 



For the ideal Davies-Cotton this takes the explicit form 



(X,Y,Z) DC 



x + 



y + 



f 



y f 2 v„ \ f 

'(A.6) 

It is straightforward to numerically calculate the image cen- 
troid Q, rj) and the resolution (A£, Arj) of such a telescope and 
to estimate the contribution of various terms to the resolution 
with a Taylor development of X and Y in terms of x, y and (f>. 
The development of terms of the form x i y J ^> k , with i + j < 5 and 
k < 3 was found to be sufficient for a percent precision in the 
resolution parameters. 

The tangential and sagittal barycenter in the focal plane (the 
image centroid) for a uniform beam on the primary surface are 
given by 



r 1/2 r 

J r dr J d6 Y(x = r cos 8, y = r sin 9) 



€ = 



V = 



j: l2 r d rj d 



rl/2 r 

J r dr J d# X(x = r cos 6, y = r sin 0) 



j: l2 r d rj d 
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Table A.l: Table of coefficients c l j for Eq. |A.9| and |A.10| 

and the corresponding resolution in term of root-mean-square: 
J^ 1/2 rdr Jd0(Y(x = rcos6»,y = rsinf?) ~lf 



An 2 = 



£ /2 r d rj d e 



J^ 1/2 rdr j d0(X(x = rcos<9,y = rsinf?) - rj) 2 



£' 2 rdrjde 



(A.8) 



The upper integration bound in r originates from the fact that 
the optical system was scaled to meet a reflector mirror of d = 1, 
hence r= 1/2. 

Taylor development of the above formulas brings the desired 
result 



at 



Arf 



U ■ 
1 c 1 
24 Zj f j Y 



(A.9) 



(A. 10) 



Coefficients given in Table A. 1 For the ideal Davies-Cotton, 

c ( . . c?. - 

only leading terms jj-cf)' and jj-tj)' are retained, at / = 1 and 
a maximum off-axis angle of the incoming rays of (p max =5°, 
ij terms such that cf .07(2/-' A£) > 1(T 3 and c'j .cp' /(2fJ Atj) > 

10 : . 

As we apply the same conditions to spherical prime-focus 
design, less terms are present at higher order, i.e. less spherical 
aberration. To mirror the result for the ideal Davies-Cotton, 
several non leading terms are added giving a consistent picture 



(A. 7) for both developments as shown in Table A. 1 
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Comparisons. A comparison between exact (numerically cal- 
culated) results and the presented limited Taylor development 



for both designs is presented in Fig. A. 7 
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Figure A. 7: Comparison between the exact (solid), i.e. numerically calculated, 
resolution parameters Af and At; and the result of the limited Taylor develop- 
ment (dashed). The inset shows the ratio of both for / = jl.0, 1.4, 2.0} (solid, 
dashd, dotted). 



In lfl6l and fP71 . 3 rd order developments for the Davies- 
Cotton and sherical mirror, respectively, have been discussed. 
The obtained coefficients are repeated here for completeness in 
table |A.l| Both solutions show up to 20% fractional error dArj, 
e.g. at / = 1. 

While in ifTTl . terms x'yj(f> k in the development were kept 
only to the 3 rd order, i.e. i + j < 3, here terms were kept up to 
i + j < 5 and k < 3. 

At the expense of the introduction of more terms, conse- 
quently, the precision of the presented development is about ten 



times better as illustrated in Fig. A. 8 



Obscuration. In the above considerations, the shadow of a pos- 
sible detector in the focal plane has been neglected. By chang- 



ing the lower bound from r = to r = 5°-/ in expressions A. 8 
obscuration can easily be quantified. Fig. A. 9 shows that ob- 



scuration degrades the resolution parameters by about 1.5% at 
/ = 1 up to about 6% at / = 2, 5° off-axis. 





/=1.4 























> [degl 



Figure A. 8: Comparison between the exact, i.e. numerically calculated, reso- 
lution parameters Af and A^ and the result of the limited Taylor development 
from fr7landfT6l. 



Appendix A.2. Parameterization for a tessellated reflector 

A realistic implementation of the non-constructable Davies- 
Cotton telescope consists in introducing a reflector made of 
multiple individual spherical mirrors. The tessellation number 
is defined as the the number N of mirrors in the diagonal. In the 
limit N = oo, it is identical to the ideal Davies-Cotton design. In 
practice, is a number <30. 

Simulation. The effective parametrization is presented as a cor- 
rection to the limited Taylor development derived earlier. The 
correction is implemented through ray-tracing simulation per- 
formed with the MARS software (described in [18] Oil ED), 
which do fully reproduce the results obtained earlier in the case 
of a spherical and ideal Davies-Cotton reflectors. Although an 
ideal Davies-Cotton reflector cannot be build in reality, it can be 
simulated easily. Simulations enable the use of arbitrary tessel- 
lation, since analytical solution being not very well suited for 
this task. 

For the simulation the following properties have been used: 

• Individual mirrors are hexagonal. For symmetry reasons 
each hexagon is rotated by 15° against the x-/y-axis 

• The mirrors are fixed on a hexagonal grid in the jc/y-plane 
with spacing d 

• The diameter of the individual mirrors is defined as nd 2 = 
nD 2 / V3, with n being the total number of mirrors in the 
system 

• Their center is located on a sphere around the focal point 
(this corresponds to foe for the ideal Davies-Cotton) 

• The focal length F of each mirror is equal to the radius of 
the sphere, and therefore equal to the focal length of the 
system 

• The overall shape of the reflector is also hexagon like 

• The tessellation number jV is the number of mirrors in the 
diagonal 

• Each mirror is oriented to a virtual point in 2F (this corre- 
sponds to «dc for the ideal Davies-Cotton) 
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Figure A. 9: Resolution parameters Af and Arj of the exact results, i.e. numeri- 
cally calculated, with (dashed) and without (solid) obscuration. Inset shows the 
ratio of both for / = { 1 , 1 .4,2} (solid, dashed, dotted). 



The small effect of obscuration by the focal instrumenta- 
tion is neglected 

The mirror surface is assumed to be ideal 



An example for such a reflector is given in Fig. A. 10 



Simulations have been carried out for N between 1 and 79 in 
steps of two, in the range 1 < F/D < 2 in steps of 0. 1 and for 
rays with off-axis angles comprised in 0°< <p <6.5° in steps of 
0.5°. 

Empirically, it could be found that introducing a dependence 



on the tessellation number, the formulas given in Appendix 



A. 1 for the spherical mirror and the ideal Davies-Cotton mir- 
ror could be unified. For this, a linear dependence at th order 
in (f>, and a quadratically at order <p' for i + 0, has been intro- 
duced. Additionally, an effective rescaling f e s = f/w is needed 
to reach an accuracy about 5% in the whole simulated range. 
The root-mean-square of a tessellated Davies-Cotton can then 
be written as 
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Figure A. 10: An example of a reflector layout for N=9 mirrors on the diagonal. 
This also includes N={3, 5, 7} removing the rows of outer mirrors consecu- 
tively. 



The coefficients s,j and dij are the ones given in Table A. 1 for 
the spherical and the ideal Davies-Cotton, respectively. 

Average correction factor for F/D=[1 ..2] 
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Figure A. 11: The scale factor w as defined by formula |ATT| determined from 
fits to the simulated point-spread function. Each point is the average of all cor- 
rection factors obtained for the simulated range of F/D. The error bar denotes 
the spread. The two curves are the tangential (upper curve, solid) and sagittal 
(lower curve, dashed), respectively. 

The tessellation number N can here be interpreted as a pa- 
rameter describing the transition from a single spherical mir- 
ror to an ideal Davies-Cotton reflector. The rescaling factor w 
can be interpreted as the deviation of the shape from the ideal 
case. Its value was determined by minimizing the residual, i.e. 
X 2 , between simulated point-spread function and approximated 
root-mean-square. The differences of the sagittal and tangential 
residual are minimized independently for each N. Its value is 



depicted in Fig. A. 11 The introduction of this scale factor ef- 
fectively reduced the residual from a maximum of 12% to less 
than 5% for tessellation numbers smaller than 40. Fig. A. 12 



(A.ll) 



shows the distribution of the residuals for different tessellation 
numbers. 

Note that for the case N — 1 the simulated single mirror 
is of hexagonal shape while the analytical approximation de- 
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Figure A. 12: Relative residual between the result of Eqs. | A. 1 T] and the simulated point-spread function after application of the correction factor w in the range 
F/D = [1,2] and o=[l,6.5]. The black lines corresponds to the median of the distribution. The gray shaded areas to 68%, 95%, and 100% of the distribution. The 
higher deviations at high tessellation ratios are located around small field-of-views (a < 3.5°) and small F/D < 1.3. The left and the right plot show the sagittal and 
the tangential component, respectively. 



scribes a disc-like mirror. For cases N » 1 the properties of the 
simulated reflector converge to the ideal Davies-Cotton. While 
the presented development was calculated for a disc shaped re- 
flector, here, the simulated Davies-Cotten converge to an ideal 
hexagon. Consequently, in both cases the rescaling factor is 
expected to be different from unity. 

In general, it is not expected to obtain a perfect match be- 
tween the analytical approximation and the simulation, because 
simulations will always take into account effect which cannot 
be easily described analytically, like rays lost between individ- 
ual mirrors. 



References 



From Eqs. A. 1 1 it is evident that for rays with small incident 
angles the point-spread function is dominated by the th -order 
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